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ABSTRACT 

Strong-Field Electrodynamics (SFE) is Maxwell theory with a certain Lorentz- 
covariant Ohm's law which uses only the electromagnetic degrees of freedom. We show 
that SFE is semi-dissipative: while the dissipation rate of the electromagnetic energy is 
non-negative, it can be exactly zero for non-trivial electromagnetic fields. 

It appears that SFE is well-defined for arbitrary electromagnetic fields. It should 
be possible to calculate the dissipative pulsar magnetosphere and resolve the magnetic 
separatrix using SFE. 

We show that SFE reduces to Force-Free Electrodynamics (FFE) in the large con- 
ductivity limit. In the regions where the ideal FFE 4-current is space-like, SFE predicts 
small dissipative corrections. In the regions where the ideal FFE 4-current is time-like, 
SFE predicts a zero correction. This indicates that bright pulsars radiate primarily 
from the magnetic separatrix. 



1. Introduction 

It has been argued that large-scale electromagnetic fields of bright Crab-like pulsars are prop- 
erly described by Force- Free Electrodynamics (FFE) (Gruzinov 2007 a). This allows to calculate the 
ideal pulsar magnetosphere (Contopoulos, Kazanas, Fendt 1999, Gruzinov 2005, 2006, Spitkovsky 
2006). 

FFE is dissipationless. In FFE the entire spin-down power is carried by the large-scale Poynting 
flux. Formally, FFE predicts zero radiation. This appears to be a good approximation for Crab- 
like pulsars - observed luminosity is indeed much smaller than the estimated spin-down power. 
Radiation should come as a small dissipative correction to FFE. 

But of course, we must find a way to calculate these small dissipative corrections, because most 
of the pulsar data concerns radiation. One way to approach this problem is to first add dissipation to 
FFE and calculate the dissipative pulsar magnetosphere (and then find the microscopic mechanism 
of the finite conductivity and the radiation counterpart of the dissipation) . 

Strong-Field Electrodynamics (SFE), a dissipative extension of FFE, has already been formu- 
lated (Gruzinov 2007 b). It was shown that in the high conductivity limit, in the regions where the 
ideal FFE 4-current is space-like, SFE reduces to FFE with small dissipative corrections. 
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On the other hand, in the regions where the current is time-like we found finite-time singular- 
ities. Because of the singularities, we erroneously concluded that SFE reduces to FFE only in the 
space-like current regions, and fails in the time-like current regions. 

The presence of singularities and the failure of the model are not equivalent 0. The singularities 
can be physical, like shocks and tangential discontinuities in hydrodynamics and current layers in 
FFE. Indeed we will argue in this paper that SFE (regularized by an arbitrarily small diffusion) is 
well-defined everywhere. One may directly use SFE to study dissipation in FFE. In particular, it 
should be possible to use SFE to calculate the dissipative pulsar magnetosphere. 

It remains to be seen if SFE will be useful in calculating the pulsar radiation, but one thing is 
already clear. If SFE gives a correct description of dissipative corrections to FFE, the emission of 
bright pulsars must primarily come from the magnetic separatrix, where the 4-current is space-like 
(and infinite in the FFE approximation). 



2. FFE and SFE equations 

FFE describes electromagnetic fields of special geometry, with electric field smaller than and 
perpendicular to the magnetic field. In 3+1 formulation, FFE is defined by Maxwell equations with 
the ideal Ohm's law 

„ „ n „ „ n . . (B-VxB-E-VxE)B + (V-E)ExB 
8 t B = -VxE, 9tE = VxB-j, j = ^ ^ i '- . (1) 

Covariant form of FFE is 

d u F^ = -f, F^j u = 0. (2) 
F is the electromagnetic field tensor. 

FFE equations come from the following physics. At each event go to a good frame, where E 
is parallel to B. In FFE, in a good frame, the electric field actually vanishes and the current flows 
along the magnetic field. This physical formulation of FFE has been used to argue that FFE works 
for Crab-like pulsars. 

SFE describes electromagnetic fields of arbitrary geometry. In 3+1 formulation, SFE is defined 
by Maxwell equations with the Ohm's law 

. _ pE X B + (p 2 + 7 2 ( r%2 ) l/2 (jjoB + EpE) 

J " B^ + El ' [6) 



where 



B 2 -E 2 =B 2 -E 2 , B E = E B, E >0, 7 2 = (4) 

+ & 



1 Here I have to quote Jeremy Goodman, "...I'm not sure that this constitutes a "failure" of SFE. ...there are 
well-defined weak solutions, i.e. solutions with discontinuities 
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p = V • E is the charge density, and the conductivity scalar a = o-(Bq, Eq) is an arbitrary function 
of the field invariants. The covariant formulation of SFE is 

d v F»» = -f, BoF^ju = EoF^j v , j^f = -a 2 E 2 . (5) 

F is the dual tensor. 

SFE equations come from the following physics. At each event, the best of the good frames 
exists, where E is parallel to B and the charge density vanishes. In the best frame, the current aE 
flows along the magnetic field. 



3. SFE in action 

3.1. Is SFE well-defined? 

The SFE Ohm's law is dissipative, E • j > - the electromagnetic energy damping is positive 
or zero in any frame. One can also show that ultraviolet modes propagate on or within the light 
cone. These good properties should be taken as an indication that SFE is well-defined. 

However, in (Gruzinov 2007b) we found finite-time singularities for certain initial fields (infinite 
second derivative of the electric field, to be discussed in §3.3). We erroneously concluded that SFE 
fails for such initial fields. 

Here we show that these singularities are physical, like shocks in compressible hydrodynamics. 
It appears that SFE is actually well-defined. We have confirmed this by a large number of numerical 
simulations of SFE, which showed no pathologies, in the sense that the evolution of the field was 
independent of the regularization. 



3.2. Regularized SFE 

The most straightforward regularization was used in Maxwell equations (realized as the Lax 
scheme in numerical simulations): 

t B = -V x E + 0- AB, d t E = V x B - j + • AE. (6) 

The SFE Ohm's law needs regularization only in the regions where both invariants of the electro- 
magnetic field vanish. We have used the following regularization when calculating Bq and Eq: 



B 2 -E 2 + J (B 2 — E 2 ) 2 + 4(E • B) 2 
Bl = > ^ } - + 0, (7) 

and then calculated Eq and Bq from unmodified expressions 

E = ^B 2 -B 2 + E 2 , Bq = sign(E ■ B)^5g. (8) 

The regularization of the sign function is arbitrary. 
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Fig. 1. — Charge relaxation in ID SFE. The electric field is shown at the initial moment (sinusoidal) 
and at three later moments (lower amplitude curves). Singularity immediately forms at zero electric 
field. Singularities also form at the extrema of the electric field. 



- 5 - 



3.3. Charge relaxation in ID 

To show both regularizations in action, consider one-dimensional charge relaxation. That is 
take an initial field B = 0, E = (E(x),0, 0). Then only the x-component of the electric field E will 
be non-zero, and SFE equations reduce to 

E = -sign(E)\J E' 2 + a 2 E 2 , (9) 

where dot and prime are the time and x derivatives. Here and below, we take a = const. 

It is easy to show that equation ([9]) gives finite-time singularities near the extrema of E. Say 
there is a positive maximum of the initial field at x = 0. Then at all time t > 0, the Taylor 
expansion near x = is E = a(t) — b{t)x 2 + o{x 2 ) with positive a, 6, and equation ([9]) gives 

b 2 

a = -aa, 6 = 2 ab. (10) 

aa 

It follows that b becomes infinite in a finite time. At this moment the second derivative of E 
becomes infinite. 

But so what? The regularized (according to §3.2) SFE equations keep working as shown in fig. 
1. In fact, we also see that singularities immediately form near the electric field zeros. The electric 
field zeros develop into zero-field regions bounded by electric field singularities (again infinite second 
derivative of the field). The electric field dies out in a finite time, even if the conductivity is zero. 
In SFE, the damping rate 

q^j-B = E ^p 2 + 1 2 a 2 E 2 (11) 
can be positive even for zero conductivity a. 



3.4. ID SFE. The FFE limit 

Generic ID SFE describes electromagnetic fields B = (0, B y ,B z ), E = {E x ,E y ,E z ), where all 
components are functions of x and t. We ran many simulations of ID SFE. It appears that SFE 
(regularized according to §3.2) works in ID. 

It is interesting to see what happens to ID FFE equilibria in the SFE description. The field 

E y = 0, E z = 0, B 2 y +B 2 z -E 2 x = const > (12) 

is an FFE equilibrium. We used the field (|12p as an initial condition for numerical SFE. It was 
found that equilibria with the space-like 4-current 

j 2 = E' 2 -B' 2 -B' 2 <0 (13) 

decay, and the decay rate decreases with the increasing conductivity a. It was also found that 
equilibria with the time-like 4-current, j 2 > 0, are not damped at all. 
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The origin of this semi-dissipative behavior is clear. In the space-like case, a non-zero parallel 
electric field is needed to sustain the current. The dissipation rate q is then positive. The dissipation 
rate decreases for increasing a, because for a given space- like current, Eq = \/—j 2 /o- decreases. 

In the time-like current regions, the parallel electric field oscillates around zero, giving a time- 
like time-averaged current. The amplitude of these oscillations is infinitesimal, and the dissipation 
q vanishes. 

3.5. 3D SFE 

We ran many simulations of full 3D SFE in a periodic cube, with resolution up to 100 3 . It 
appears that SFE (regularized according to §3.2) works in 3D too. 

We checked, for example, that a space-like FFE equilibrium E = 0, B = (sin y — cos z, sin z — 
cos x, sinx — cos y) is damped and the damping rate goes down for increasing a. We also observed 
formation of zero-field regions and infinite second derivatives for a spherically symmetrical pure 
electric field. 

I thank Jeremy Goodman who pointed out that SFE works better than I originally thought. 
This work was supported by the David and Lucile Packard foundation. 
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